
program ttrajd

! testovani prevodu Julinaskeho data

use trajd

real(double) :: y,m,d,jd

y = -605
m = 10
d = 11.0
jd = datjd(y,m,d)
write(*,*) y,m,d,jd
call jdat(jd,y,m,d)
write(*,*) jd,y,m,d

end
